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ABSTRACT 



1 INTRODUCTION 



The large scale structure of the Universe, as depicted from 
galaxy surveys, weak lensing mapping and numerical simu- 
lations, shows a web like three dimensional structure. There 
are three features that can be generally observed. First, most 
of the volume resides in underdense regions; second, most of 
the volume is permeated by filaments; third, the densest 



A dynamical classification of the cosmic web is proposed. The large scale en- 
vironment is classified into four web types: voids, sheets, filaments and knots. The 
classification is based on the evaluation of the deformation tensor, i.e. the Hessian 
of the gravitational potential, on a grid. The classification is based on counting the 
number of eigenvalues above a certain threshold, A t h, at each grid point, where the 
case of zero, one, two or three such eigenvalues corresponds to void, sheet, filament or 
a knot grid point. The collection of neighboring grid points, friends-of-friends, of the 
same web attribute constitutes voids, sheets, filaments and knots as web objects. 

A simple dynamical consideration suggests that A t h should be approximately unity, 
upon an appropriate scaling of the deformation tensor. The algorithm has been applied 
and tested against a suite of (dark matter only) cosmological N-body simulations. In 
particular, the dependence of the volume and mass filling fractions on Ath and on the 
resolution has been calculated for the four web types. Also, the percolation properties 
of voids and filaments have been studied. 

Our main findings are: (a) Already at Ath = 0.1 the resulting web classification 
reproduces the visual impression of the cosmic web. (b) Between 0.2 < Ath S 0-4, a 
system of percolated voids coexists with a net of interconected filaments. This sug- 
gests a reasonable choice for Ath as the parameter that defines the cosmic web. (c) The 
dynamical nature of the suggested classification provides a robust framework for incor- 
porating environmental information into galaxy formation models, and in particular 
the semi-analytical ones. 

| |White fc' R.ces 1978). It follows that the study of such en- 
vironmental dependence should commence with the effort 
to understand the formation of DM halos in the context of 



the cosmic web dGao et al.||2005| lAvila- Reese et al. 1120051 



Maulbetsch et al. 2007b. This motivates us to search for a 



clumps arc located at the intersection of filaments (Bond 
et al.|1 996). This motivates a classification of the cosmic web 



into at least three categories: voids (underdense regions), fil- 
aments and knots (densest clumps). 

There are clear evidences for the correlations of the ob- 
served properties of galaxies with the environment. We have 
for instance the morphology density correlation, where el- 
liptical galaxies are found preferentially in crowded environ- 



ments, and spiral galaxies are found in the field (Dressier 



robust and meaningful method to classify the different en- 
vironments in numerical simulations. Such a classification 
should provide the framework for studying the environmen- 
tal dependence of galaxy formation. 

Translating the visual impression into an algorithm that 
classifies the local geometry into different environments is 
not a trivial task. A somewhat less challenging, yet very 
closely related, task is that of identifying just the voids out 
of the cosmic web. A thorough review and comparison of dif- 
ferent algorithms of void finders has been recently presented 
in Colberg et al. (20081. The void finders can be classified 



19801. The same kind of correlation can be found in terms 
of the colors of the galaxies ( |Blanton et al.|2005| |. 

According to the current paradigm of structure forma- 
tion galaxies form and evolve in dark matter (DM) halos 



according to the method employed. Most are based on the 
point distribution of galaxies or dark matter (DM) halos and 
some on the smoothed density or potential fields. Some of 
the finders are based on spherical filters while others assume 



no inherent symmetry (Brunino et al. 2007 Colberg et al 



Forero-Romero et al. 



2005) jGottlober et al.|[2003lJArbabi-Bidgoli fc Miiller|2002 



Neyrinck|2008||Platen et al.|2007||Plionis fc Basilakos|2002 



Shandarin et al.|2006 1 

It should be emphasized that an environment finder 
should be evaluated by its merits and it cannot be labeled as 
correct or wrong. A good algorithm should provide a quanti- 
tative classification which agrees with the visual impression 
and it should be based on a robust and well defined numer- 
ical scheme. Yet, it is desirable for an algorithm to be based 
on an analytical prescription so that its outcome can be cal- 
culated analytically. Also, simplicity is always very highly 
desired. 

A variety of approaches have been employed in the clas- 
sification of the cosmic environment into its basic elements. 
The simplest way is based on the association of the envi- 
ronment with the local density, evaluated with a top-hat 
filter, say, of some width ( Lemson fc Kauffmann|1999 1. The 
density field can be analyzed in a much more sophisticated 
and elaborated way. This is the case of the web classification 
based on the multi-scale analysis of the Hessian matrix of 
the density field ( |Arag6n-Calvo et al.||2007 1 or the skeleton 



analysis of the density field (Novikov et al. 2006 Sousbie 
et al.| [2008). Both methods classify the cosmic web by pure 



geometrical tools applied to the density field. A very differ- 
ent approach is done within a dynamical framework in which 
the analysis of the gravitational potential is used to classify 
the web. This has been inspired by the seminal work of |Zel-| 
dovich (19701 that led to the "Russian school of structure 



formation" (e.g. Arnold et al. (19821; Klypin & Shandarin 



( 1983 1 ) . The quasi-linear theory of the Zeldovich approxi- 



mation predicts the existence of an infinitely connected web 
of pancakes (i.e. sheets), filaments and knots. This morpho- 
logical classification is based on the study of the eigenvalues 
of the deformation tensor, namely the Hessian matrix of the 
linear gravitational field. 

A recent application of the Zeldovich-based classifica- 



tion has been provided by Lee & Lee ( 2008 1 who used a 



Wiener filter linear reconstruction of the local density field 
and evaluated the linear deformation, and hence also the 
shear, tensor on a grid. The cosmic web has been classified 
according the structure of the shear tensor. 

A different approach has been followed by H ahn et al.| 



(20071 who suggested that the full non-linear gravitational 



potential should be used for the geometrical classification. 
Apart from the difference between the linear and the non- 



linear potential, both Hahn et al. (20071 and Lee & Lee 
( |2008[ ) are using the same classification. Namely, the Hessian 
of the gravitational potential is evaluated on a grid and its 
eigenvalues are examined locally. A grid point is classified 
as a void, sheet, filament or knot point if the number of 
eigenvalues greater that a null threshold is zero, one, two or 
three. 

The algorithm presented here is an extension and im- 



provement of the one suggested by Hahn et al. (20071. The 



extension is represented in the selection of a new free param- 
eter with a dynamical interpretation. As such it provides 
a classification of local environment. Namely each spatial 
point is flagged as belonging to a either a void, sheet, fila- 
ment or a knot point. It is the collective classification of all 
points in space which gives rise to the geometrical construc- 
tion we call the cosmic web. This opens the door for defin- 
ing voids, sheets, filaments and knots as individual objects. 



Each object is defined as a collecion of connected points 
having the same environmental attribute. We use a friends- 
of-friends (FoF) algorithm to find connected sets of points in 
simulations. Having defined the web objects the statistical 
and dynamical properties of these can be readily studied. 
Here we shall focus on analyzing the statistical properties of 
the voids and filaments, aiming to better constrain the value 
for the free parameter we introduce. We perform as well a 
statistical study to asses effect of the cosmic variance on our 
conclusions. 

This paper is organized as follows. The web classifica- 
tion scheme is described in Sj2] The N-body simulation used 
in the paper and the numerical implementation of the web 
classification are presented in ^3] ^describes the main prop- 
erties of the cosmic web and in particular its dependence on 
the smoothing scale and the free parameter of our classifica- 
tion scheme. ^5] concentrates on the properties of the voids 
sector of the cosmic web. Sj6] studies the fragmentation of 
filaments in order to give a confidence interval to the free 
parameter that was introduced. In ^7] we revisit the sections 
Sj4]and [j5]to study the effect of cosmic variance. The paper 
concludes with a general discussion and a summary of the 
main results of the paper (§81. 



2 WEB CLASSIFICATION 



Hahn et al. (20071 have recently suggested a new dynami- 



cal classification of the cosmic web. The basic idea of their 
approach is that the eigenvalues of the deformation tensor 
determine the geometrical nature of each point in space. 

The deformation tensor, T a p, is defined by the Hessian 
of the gravitational potential <fr. 



T a fi 



d 2 <t> 



(1) 



The definition of the deformation tensor Q explicitly as- 
sumes that the matter density field is known and that it is 
smoothed with a finite kernel, or otherwise the derivatives 
are not defined. For simplicity the smoothed density field is 
defined over a (Cartesian) grid. 



Hahn et al. (20071 considered the three eigenvalues of 
the deformation tensor, Ai > A2 > A3 , and classified a grid 
point according the number of positive eigenvalues at that 
point. Namely, a void point corresponds to no positive eigen- 
values, a sheet to one, a filament to two and a knot point 
to three positive eigenvalues. The sign of a given eigenvalue 
at a given grid node determines whether the gravitational 
force at the direction of the principal direction of the corre- 
sponding eigenvector is contracting (positive eigenvalue) or 
expanding (negative). 

|Hahn et al.| ( |2007[ ) provided a very attractive approach 
to the web classification problem. It is based on the dy- 
namical nature of the web, and so it easily lends itself to 
a theoretical analysis. The ease of its application to cosmo- 
logical simulations opens the door for a new framework for 
associating the properties of galaxies and dark matter with 
environment, as defined by the web classification. 



Hahn ot al. 



1 2007) call the T a p tensor the 'tidal tensor'. Usually 



the tidal tensor is defined as the traceless part of T a p. 
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Close inspection of the Hahn et al. ( 2007 I classification 



scheme and its results reveals its shortcomings. The volume 



filling factor of Hahn et al. ( 2007 1 voids is very small. For 
their minimal smoothing scale, namely the highest resolu- 
tion, the voids occupy only 17% of the simulated volume. 
This stands in contrast to the visual impression of voids in 
the actual universe and in simulation, where voids seems to 
occupy most of the volume but contain only a small fraction 
of the galaxies (in observations) or matter (in simulations). 
Furthermore, the Hahn et al. (20071 classification does not 



reproduce the visual perception of the cosmic web. 

It is easy to understand the inability of the |Hahn et al.| 
( |2007[ > approach to reproduce the visual impression. The 
web classification is based on the algebraic sign of the eigen- 
values of the deformation tensor, namely the number of 
eigenvalues larger than a threshold value of zero. It fol- 
lows that if an eigenvalue is only infinitesimally positive, the 
scheme assumes that the local neighborhood of the given 
grid point collapses along the corresponding eigenvector. 
Yet, the collapse proceeds over the dynamical time scale 
and if the value of the eigenvector is small enough the col- 
lapse will occur, if at all, only in the distant future. Visual 
inspection would not classify the region as collapsing at the 
present time. This leads us to consider an alternative ap- 
proach, namely asking the eigenvalues to be larger than a 
positive threshold. 

Seeing that the dimensionality of the deformation ten- 
sor is [time] -2 , its eigenvalues can be associated with the 
collapse time. It follows that the threshold value should be 
roughly determined by equating the collapse time with the 
age of the universe. In Appendix [X] we rewrite the eigen- 
vectors of the deformation tensor in a dimensionless way. 
In such a presentation one expects Ath ~ 1, where A t h is 
the threshold value. For an isotropic collapse A t h can be cal- 



culated explicitly (see Eqs. (A3 1 - (A7l). However, this pro- 
vides a very rough guide as the collapse on the web is clearly 
not isotropic. Here, an empirical approach is to be used and 
Ath is to be roughly determined from the geometrical nature 
of the cosmic web it defines. 



3 N-BODY SIMULATION, NUMERICAL 
IMPLEMENTATION AND OBJECT 
DETECTION 

We use two numerical simulation. The first assumes a 
WMAP3 cosmology ( |Spergel et al.|2007| with a matter den- 
sity Sl m = 0.24, a cosmological constant Qa = 0.76, a dimen- 
sionless Hubble parameter h = 0.73, a spectral index of pri- 
mordial density perturbations n = 0.96 and a normalization 
of o-g = 0.76. A simulation of box size 160/i _1 Mpc and 1024 3 
particles has been assumed, corresponding to a particle mass 
of 3.5 x 10 8 Mq. Starting at redshift z — 30 the evolution is 
followed using the MPI version of the Adaptive Refinement 
Tree (ART) code described in |Gottlober fc Klypin| ( [20081 ) . 
The simulation used here is actually a constrained simula- 
tion of the local universe which is to be described at length at 
the forthcoming Yepes et al paper (in preparation). This is 
an updated and higher resolution version of the constrained 



simulation presented in Klypin et al. (2003). Here the sim- 



ulation is treated as just a random one and its constrained 
nature is completely ignored. 



In order to estimate cosmic variance effects we have 
used a numerical simulation of box size l/i^Gpc. The as- 
sumed cosmology for this simulation is WMAP3-like with a 
matter density f2 = 0.27, a cosmological constant Oa = 0.73, 
a dimensionless Hubble parameter h — 0.70, a spectral in- 
dex n = 0.95, a normalization as = 0.79 and 1024 3 particles 
corresponding to a particle mass of 9.8 x 10 10 M Q . The sim- 
ulation was performed using the ART code as well. 

The analysis of the simulations proceeds as follows. The 
density field of the 160/i~ Mpc simulation is calculated from 
the particle distribution on a 256 3 grid using the Cloud-In- 
Cell (CIC) scheme, it is then smoothed with a Gaussian 
kernel of width R s and from which the deformation tensor 
is calculated directly, using an FFT solver. The deformation 
tensor is then diagonalized on the grid. The web characteris- 
tic of each grid point is determined by the number of eigen- 
vectors, at that grid node, above the threshold. It should be 
realized that the classification is local by its nature, but the 
combined effect of all grid points results in the geometrical 
construction defined as the cosmic web. 

In the l/i _1 Gpc simulation we followed the same proce- 
dure, interpolating first the density on a 512 3 grid. Once the 
environment detection procedure is done we select 6 3 = 216 
non-intersecting sub-boxes of 160/i~ Mpc on a side to test 
the effect of cosmic variance. Figure [T] presents the CIC den- 
sity field and the cosmic web evaluated at the threshold val- 
ues of A th = 0.00, 0.20, 0.40, 1.00 and 2.00. The web is pre- 
sented by a grey scale corresponding to the four web types, it 
is evaluated at a Gaussian smoothing of R s = 0.625/i -1 Mpc 
in the Figure [T] The density field and the cosmic web are 
evaluated on a plane of the CIC grid. Visual inspection of 
the density field reveals a network of voids, filaments and 
dense knots. The density field is evaluated on a thin plane 
and therefore no clean distinction can be made between the 
3-dimensional filaments and sheets. 

The web denned by A t h = 0.0 consists of many small 
isolated voids that occupy only a small fraction of the total 
area of the plane. Only as Ath increases the voids get bigger 
and connected and they become the dominant geometrical 
component of the web. 

The non-zero threshold classification provides a better 
visual match to the density field than the null case. A quali- 
tative analysis and comparison is presented in ij4] The anal- 
ysis is based on two quantities: the volume occupied by each 
web type (volume filling fraction - VFF) and the fraction 
of mass contained in such a volume (mass filling fraction - 
MFF). 

A characterization of the cosmic web is obtained by 
grouping, by a Friends-of-Friends (FoF) algorithm, neigh- 
boring grid points of a given web type into individual ob- 
jects. The resulting objects are defined as voids, sheets, fil- 
aments and knots. The FoF association proceeds in the fol- 
lowing way: the centers of the cells in the grid are used as 
the position of four different kinds of particles according to 
its web type, then a standard FoF is run over particles of 
the same kind with a linking length b — 1.1 times the grid 
length. It means that only the six closest neighbors of a 
given cell are taken into account. The FoF void detection is 
done for different simulations: the 160/i~ 1 Mpc simulation, 
the l/i _1 Gpc simulation and 216 sub-boxes of 160/i _1 Mpc 
on a side extracted from the largest simulation. To be con- 
sistent in the analysis of these three kinds of simulations, 
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Figure 1. The density field and five different environmental classifications. Upper left panel: slice of width 0.625/i _1 Mpc depicting the 
density field in the simulation smoothed over a scale of _R S = 0.625fe~ 1 Mpc, the color coding is logarithmic in the density — high density 
peaks are dark. The other panels show the environment classification using different values in for the threshold, A t h = 0.0, 0.20, 0.40, 1.0, 
and 2.0. White corresponds to voids, clear gray to sheets, dark gray to filaments and black to peaks. The general impression is that the 
non-zero values of A t h below 1.0 capture better the environment seen by eye in the density plot. 
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one should make the FoF detection without considering the 
periodic boundary conditions, as it is the right boundary 
condition for the sub-volumes extracted from the l/i _1 Gpc 
simulation. 

The FoF algorithm is used to detect voids for differ- 
ent threshold values in the different simulations which are 
smoothed with the same physical scale R B = 1.95/i _1 Mpc 
in order to allow a fair comparison between them. In Sj5] a 
detailed analysis on these voids is presented, paying special 
attention to its percolation properties as the threshold rises. 
The same kind of analysis is performed for the filaments in 
fjfi] but only on the l/i _1 Gpc simulation. 



4 VOLUME AND MASS FILLING FRACTIONS 

The web classification depends on two parameters that de- 
termine the environment. The first is the smoothing scale _R S , 
and the second is the threshold for the eigenvalues Ath. The 
dependence of the volume and mass filling fractions (VFF 
and MFF, respectively) on these two parameters is studied 
here. This is done for the four web types. 

The VFF and MFF are measured for every kind of en- 
vironment. First by fixing A t h = and varying the smooth- 
ing scale R B between 0.625/i -1 Mpc and 12.4/i _1 Mpc. Later, 
fixing the smoothing scale to 1.95/i _1 Mpc and varying Ath 
between and I with steps of 0.1. The results of these two 
kinds of cuts are shown in Figure [2] The VFF and MFF of 
the case A t h = 0.0 and 1.0, with two different smoothing 
scales, R B = 0.625/i _1 Mpc and R B = 1.95/i _1 Mpc are pre- 
sented in Table [I] The evolution of the filling fraction with 
Ath = reproduces the asymptotic results expected for large 
smoothing scales, which is 0.42 for sheets and filaments and 
0.08 for voids and knots ( Doroshkevich|1970 ). 

The most striking feature that emerges from Figure [2] 
is that of the strong dependence of VFF and MFF of the 
voids on A t h- This is to be contrasted with the other web 
types which show quite a similar behavior with the change 
of the smoothing and the change of the threshold. It follows 
that the voids can serve as a sensitive monitor and indica- 
tor of the cosmic web. For the case of a null threshold the 
dependence on the smoothing length is weak. The increase 
of Rs corresponds to a transition to the linear regime where 
the density field is closer to be Gaussian. 

Then we perform an exploration of the space A t h-i? s by 
making the environment classification with steps of AA t h = 
0.1 for each smoothed density field. We measure the filling 
fractions for each couple i? s -A t h- The results are shown in 
Figure [3] The important result here is that in the range of 
explored values, the threshold on the eigenvalues is more 
important that the smoothing scale to fix the volume filling 
fraction. 

Inspection of the behavior of the VFF and MFF in the 
R B - A t h plane reveals that voids on the one hand and the 
sheets and filaments on the other hand are complementing 
one another. Meaning that the VFF and MFF for voids in- 
crease with both R B and Ath at expenses of the volume and 
mass in sheets and filaments. 

This provides another evidence for the distinct nature 
of voids. The VFF of the voids, sheets and filaments clearly 
show the significance of the Ath J; 0.2 case. For that thresh- 
old the contours lines are almost vertical, implying that the 



VFF of these web elements are almost independent of R B . 
The MFF, on the other hand, distinguishes very easily dif- 
ferent values of the smoothing scale Rs. 



5 PERCOLATION OF VOIDS 

Given that each grid node is flagged as belonging to a void, 
sheet, filament or a knot one can define new objects made 
by connecting grid nodes of the same type. In particular 
we proceed here to define a void as an object made of the 
collection of neighboring void type grid nodes. Neighboring 
nodes are associated by the FoF algorithm to form individual 
objects as described at the end of [J3] 

Here we focus on the analysis of the statistical distribu- 
tion of the sizes of the voids and their percolation. The em- 
phasize on voids does not stem only from the extensive work 
done on their properties ( Colberg et al.p008 l but also be- 
cause they constitute the most sensitive gauge of the cosmic 
web and its dependence on the threshold of the eigenvalues. 
Filaments have received attention as well in the literature 
( Novikov et al.[ 2006), in the next section we present a brief 
analysis of their properties. 

A simplified characterization was already performed in 
!^4]by measuring the volume and mass occupied in the void 
environment. Here, the dependence of the number of voids 
and their percolation properties on A t h is examined. The 
percolation is quantified now by the fraction of volume of 
the most voluminous void to the total volume occupied by 
all voids. 

The void identification is performed here at a fixed 
smoothing scale R B — 1.95/i _1 Mpc and by varying the eigen- 
values threshold in the range of < A t h < 0.3. In the Figure 
[I] (left panel) one can see that the number of voids in the 
simulation roughly decays exponentially as a function of the 
threshold parameter A t h, N void = No exp(— A t h/Au), where 
Nq is the number of voids for the null threshold and A_d. 

Figure [4] (right panel) presents the fraction of the vol- 
ume of the largest (in volume) void to the total volume of 
all the voids. A transition occurs between 0.1 < Ath S; 0.2 
where it jumps from a ratio < 0.1 to > 0.9. Note that the 
percolation starts at the stage in which the VFF of voids 
is only ~ 25%. It follows that in spite of the small VFF 
obtained for Ath ~ 0.1 the voids start to coalesce and form 
one super-void which encompasses 90% of the total volume 
of voids when the void VFF reaches 60%. As we will show 
in SjT] this transitional scale is dependent on the particular 
simulation under consideration. 



6 FRAGMENTATION OF FILAMENTS 

Along with voids, filaments have received some attention in 
the literature as a natural way to probe large scale structure 
( jSousbie et al.|2 008). A study of the filaments is performed 
here, focusing on their percolation dynamics. 

In the opposite sense to voids, filaments start to frag- 
ment as the threshold A t h is raised. This is clear in Figure [I] 
which shows that the filaments dissappear as A t h increases, 
making room for the growing voids. 

Consequently, one could expect that for some range of 
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5 10 0.0 0.5 1.0 

R s [h- 1 Mpc ] A, th 

Figure 2. Upper panel: the volume filling fraction of peaks, filaments, sheets and voids as a function of the smoothing scale R s for 
Ath = 0.0 (left), and as a function of A t h for R s = 1.95/i _1 Mpc (right). Continuous line: voids, dashed: sheets, dotted: filaments, 
dotted-dashed: knots. Lower panel: same as the upper panel but for the mass filling fraction. 
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Web type 


Volume 
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Mass 


Volume 


Mass 


Volume 


Mass 


void 


0.16 


0.02 


0.76 


0.28 


0.13 


0.03 


0.82 


0.47 


sheet 


0.60 


0.27 


0.18 


0.25 


0.56 


0.32 


0.14 


0.25 


filament 


0.24 


0.54 


0.05 


0.35 


0.28 


0.52 


0.04 


0.22 


knot 


0.01 


0.16 


5.0e-3 


0.12 


0.01 


0.11 


2.8e-3 


0.06 



Table 1. Volume and mass filling fractions for the four web types for two different smoothing scales, R s = 0.625/i -1 Mpc and R s = 
1.95/i _1 Mpc. For each smoothing two extreme values of the threshold are used, A t h = 0.0 and X t h = 1-0. The volume filling fractions are 
similar for the same values of the threshold A t h regardless of which smoothing scale is used. 
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Figure 3. Upper row. Isocontours for the volume filling fractions for the different kinds of environment, in the plane R s -X t ^. Lower row. 
Isocontours for the mass filling fractions for the different kinds of environment, in the plane R B -\ t ^. In both rows, from left to right: 
filaments, sheets and voids. The label on each isocontour is the percentage of the total volume filled by a given kind of environment. 




Figure 4. Percolation results. Left panel: number of voids as a function of the threshold in eigenvalues normalized to the number of 
voids at A t h = 0.0. A fixed smoothing of R B = 1.95/i Mpc is used here for all the simulations. The thick line shows the results for the 
lh~ x Gpc simulation, the thin line shows the results for the 160/i -1 Mpc simulation. The gray lines show the results for the sub-volumes 
extracted from the lh~ x Gpc simulation. Right panel: fraction of the total void volume occupied by the most voluminous void. The line 
coding is the same as in the left panel. Results from the large l/i —1 Gpc differ greatly from the results in the simulation 160/i -1 Mpc, 
nevertheless the later is consistent within the scatter deduced from the sub- volumes. Even when it is clear that the detailed shape of the 
curve depends on the simulation size, it seems to be a robust feature that the largest change in the super-void size is presented around 
0.1 < A th < 0.2. 



8 Forero-Romero et al. 




Figure 5. The continuous line shows the fraction of total volume 
in voids occupied by the largest void. The dashed line shows the 
fraction of the total volume in filaments occupied by the largest 
filament. Both curves refer to the lh _1 Gpc simulation. As the 
eigenvalue threshold rises the voids percolate, while the network of 
filaments fragments. The two fractions are the same at A t h = 0.25. 



Ath there are two coexisting environments: a complex of per- 
colating voids and a network of interconnected filaments, 
something close to the visual impression of the comic web. 
The percolation analysis of the filament and void networks 
can help to define an interval of Ath values where the envi- 
ronment studies would be feasible. 

We use the lft _1 Gpc simulation following the approach 
of the previous section, Figure [5] shows the ratio of the vol- 
ume of largest filament to the total volume occupied by fila- 
ments, overplotted is the same fraction for voids (presented 
in the right panel of Figure [4J. The filament fraction evolves 
from a value close to ~ 1 for A t h = 0.0 down to values close 
to zero for A t h ~ 1.5. It confirms the visual intuition we had 
of a fully interconnected network of filaments that fragments 
as the eigenvalue threshold increases. 

The percolation/fragmentation curves intersect at 
Ath = 0.25, when the two fractional volumes are ~ 97%. 
Heuristically, one can assume a given network to exist, 
namely percolate, when its fractional volume exceeds 95%. 
The web is then defined at the threshold level at which the 
voids and filaments networks coexist. This implies a thresh- 
old interval of 0.20 < A t h < 0.40. Such a heuristic approach 
is in a good agreement and matches the visual impression of 
the LSS (Figure [T]). This approximate interval should hold 
for lower smoothing scales, as the volume filling fraction (our 
gauge for the percolation dynamics) is almost independent 
of the smoothing scale for the range of thresholds consid- 
ered, as seen in Figure [3] and suggested as well in the values 
of Tabled] 



7 SIMULATION BOX SIZE AND COSMIC 
VARIANCE 

In order to quantify the influence of the simulation box size 
and the impact of cosmic variance we measure for three 



different kind of simulations the MFF and VFF for voids 
and filaments, and the percolation of voids. The first is 
the l/i _1 Gpc simulation with the density field interpolated 
over a 512 3 grid and smoothed with R B = 1.95 /i~ 1 Mpc, 
used to quantify the effect of simulation size. The second 
is the 160fe _1 Mpc simulation interpolated over a 256 3 grid 
smoothed over the same physical scale. The third kind is the 
lfo _1 Gpc simulation split in smaller non-intersecting cubes 
of 160fo _1 Mpc, used to quantify the effect of cosmic variance. 

Figure|6]shows the result for the VFF and MFF. The re- 
sults for the 160/i -1 Mpc simulation are well within the vari- 
ance calculated from the sub- volumes in the l/i _1 Gpc simu- 
lation. In the case of the MFF, the result for the 160/i _1 Mpc 
are located far from the mean value in the l/i _1 Gpc, never- 
theless it is located with in the dispersion defined by the sub- 
volumes. In general, the VFF and MFF are consistent in all 
the three kinds of simulations. The agreement is less impres- 
sive than with the VFF, perhaps due as well to the different 
values of as used in the simulations, the large simulation has 
erg = 0.79, while the small simulation has ag = 0.75 

The growth of the largest void in the simulation, Fig- 
ure [4j is very different between the three kinds of simula- 
tions. The fraction of void volume occupied by the super- 
void is very dependent on the simulation size. For the 
large l/i _1 Gpc simulation, the initial values of V ma x/Vtot 
for A t h = 0.0 are the lowest possible, this can be readily 
understood by the fact that Vtot grows with the simulation 
box size, while V max should be on the same order of magni- 
tude regardless of the simulation size. For values larger than 
Ath = 0.2 the large simulation has almost percolated into a 
single super-void, while the sub-volumes still show a large 
dispersion in their percolation behavior. The results for the 
small 160ft _1 Mpc simulation are consistent with such dis- 
persion. In spite of that, in all the three cases there is a 
clear transitional behavior starting at A t h = 0.1 and finish- 
ing around Ath = 0.2, even when the detailed behavior with 
A t h is far from being the same. 



8 CONCLUDING REMARKS 

This paper presents an improved method to identify large 
scale environment in dark matter simulations. Our scheme 
is based on the analysis of the Hessian of the gravitational 
potential generated by the dark matter distribution. The al- 
gorithm presented here constitutes an improvement on the 



scheme of Hahn et al. (20071, involving a pertinent reinter- 



pretation of the dynamics in the problem. 

The change is done through the addition of a free pa- 
rameter Ath related to the dynamical time associated with 
the collapse. This important improvement allows a more re- 
alistic treatment of the cosmic web. 

Inspection of the different plots of Figure [2] reveals the 
striking difference in the way the cosmic web, and in par- 
ticular the voids, respond to the changes in the Gaussian 
smoothing and the Ath's threshold. Keeping a null Ath and 
changing i? s we see that the non-linear evolution does not 
change the ranking of the VFF and MFF found in the deep 
linear regime (i.e., R B w 12.5/i -1 Mpc). Namely, for the null 
threshold the sheets have the highest VFF and the filaments 
the highest MFF, independent of R s . Considering the case 
of a fixed R s = 1.95/i _1 Mpc, the void VFF grows strongly 
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Figure 6. Cosmic variance effects on vofume and mass filling fractions for voids and filaments. Thick black line: simulation l/i _1 Gpc, 
thin black line: simulation 160h~ iJVlpc, gray lines: sub-volumes extracted from the x Gpc simulation. Upper panels: volume filling 
fraction of voids (left) and filaments (right) as a function of the eigenvalue threshold A t h. Lower panels: same as the upper panels for 
the mass filling fractions. 



with A t h and above Ath ~ 0.5 the voids have the highest 
VFF. The results on the VFF are extremely robust respect 
to the simulation size and cosmic variance effects. In the 
case of void MFF, it changes from the lowest one at A t h = 
to the second highest at Ath = 1.0. 

The nature of the web changes dramatically with the 
threshold and it becomes volume dominated by the voids as 
Ath increases. The MFF shows a larger cosmic variance, com- 
pared with the VFF. Also, it is equally sensitive to changes 
in the smoothing scale and the threshold value. 

The web classification provides a set of flagged points on 
a grid and the collection of neighboring grid points of a given 
environmental type, connected by a FoF algorithm, forms 
objects we call voids, sheets, etc. The statistical properties of 
the system of voids and their dependence on A t h have been 
explored here. In particular the number of isolated voids 
and their sizes have been analyzed, finding that the number 
of isolated voids roughly decreases exponentially with Ath. 



In the l/i _1 Gpc simulation the system of voids percolates 
between 0.1 < Ath < 0.2, at which the largest void jumps 
from having less than 10% to 90% of the volume occupied 
by voids. The percolation dynamics is also seen on average 
in the different smaller simulations, with a wide spread on 
the way the percolation is observed, but keeping the same 
threshold interval for the transition. 

The association of the eigenvectors of the deformation 
tensor with the collapse time, and hence with the age of the 
universe, enables in principle a theoretical determination of 
Ath- Using the spherical collapse model, calculated within 
the WMAP3 cosmology, the threshold value is A t h = 3.21. 
However, at that high value the web looks very fragmented, 
in particular the network of filaments. The application of 
the ellipsoidal collapse model (Sheth & Tormcn 2002) might 
provide a better theoretical estimate. 

Short of a "first principle" determination of A t h we re- 
sort to a heuristic approach. We look for the range of Ath over 
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which the percolated networks of the voids and filaments co- 
exist for a fixed smoothing scale of R B — 1 .95/i _1 Mpc in the 
lft _1 Gpc simulation. The voids and filaments behave in an 
opposite way in terms of the dependence of the percolation 
on the threshold value. At low A t h the voids are isolated and 
the filaments percolated and at high Ath the voids percolate 
and the filaments are fragmented. Adopting a 95% in the 
fractional volume as defining the percolation transition we 
find the web to be defined by a threshold in the range of 
0.20 < A t h 0.40. This range stands in good agreement 
with the visual impression obtained from the simulations. 

The notion of the cosmic web is not new. The filamen- 
tary structure has been extensively studied, mostly within 
the context of the Zeldovich pancakes ( Zeldovich| 1970| > . The 
role of voids has also been heavily studied and many algo- 



rithms for voids finding have been suggested (see Colberg 
|et al.| p008| ). We have been motivated by the computa- 
tional simplicity and the elegance of the Halm et al. (2007) 



approach and have modified it in a way that reproduces 
the web as it emerges from observations and simulations. 
Our main drive is to provide a simple, fast and precise tool 
for classifying the environmental properties of each point in 
space. Using the non-zero thresholding of the eigenvalues of 
the Hessian of the potential indeed provides a very efficient 
tool that can be easily applied to simulations. 

The same analysis can be performed at different red- 
shifts in the simulation, allowing the classification of envi- 
ronment as a function of time. 

The dynamical nature of the web classification implies 
that the web type might affects the dynamical evolution 
of DM halos and of galaxy formation. This might have 
profound consequences on the star formation timescale.... 
(Gayler et al. 2008). It follows that the web classification 



can be introduced into semi-analytical modeling, as a dy- 
namical tag that together with the mass of the DM halos 
dictate the the mode gas accretion onto galaxies. 

A major challenge that is still to be addressed is the ap- 
plication of the method to the distribution of galaxies. Short 
of that, the algorithm remains in the theoretical realm of 
simulations and semi-analytical modeling of galaxy forma- 
tion. 
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APPENDIX A: SETTING THE THRESHOLD 
FOR WEB CLASSIFICATION 

In this paper we use the eigenvalues of the deformation ten- 
sor, Eq.|l]), normalized in a specific way. Here we provide 
details of the normalization and give motivation for select- 
ing the threshold. 

We write the Poisson equation in the following form: 



4-rrGpS = Ai + A 2 + A :! 



(Al) 



where p is the mean matter density of the universe and S 
is the matter overdensity. One can re-scale the gravitational 
potential and the eigenvalues of the deformation tensor by 
dividing them by 4<f)Gp: 



V <h = 8 = Ai + A 2 + A 3 



(A2) 



Note that 4(f>Gp provides a natural scale to introduce 
dimensionless parameters Ai. We solve this equation numer- 
ically. 

The spherical collapse model is invoked here so as to 
get a rough estimate of A t h- The (spherical) free-fall time is 
related to the local density by 



3tt 
32Gp' 



Recalling that 



47vGp= ^ m #0, 



(A3) 



(A4) 



(where Q m is the value of the cosmological matter density 



and Ho is the Hubble constant), Eq.( Al l is rewritten as 



Ai + A 2 + A 3 = 4ivGp - AixGp. 



(A5) 



The threshold can be estimated by demanding that the 
free-fall time equals the age of the universe (to). Namely, 
the threshold separates between the principal axes that have 
collapsed by To and the ones that have not. Substituting the 
free-fall time by the age of the universe the threshold is given 
by 



3tT 

3/3(Ai)A t h = A 



-Q m Ho, 



(A6) 



where ft is a fiducial factor introduced to account for 
the deviation from local isotropy. 

In terms of the dimensionless eigenvalues the threshold 
is given by: 



Ath — 



1 

3/3(Ai) 



7T^ J_ 
4 Qrr. 



(r H y 



(A7) 



For the WMAP3 parameters used in the 160/i _1 Mpc 
simulation we have Q m = 0.24, h — 0.73 and tqHo = 0.983. 
It follows that Ath = 9.63/(3.0/3). 

Assuming the spherical collapse model as a a proxy 
to the full non-isotropic case, namely /3(\i) = 1, then the 
threshold is \ t h = 3.21. 
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